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Abstract: A first result of the search for z/„ —> v e oscillations in the OPERA experiment, 
located at the Gran Sasso Underground Laboratory, is presented. The experiment looked 
for the appearance of v e in the CNGS neutrino beam using the data collected in 2008 
and 2009. Data are compatible with the non-oscillation hypothesis in the three-flavour 
mixing model. A further analysis of the same data constrains the non-standard oscillation 
parameters 6 new and Am^ OT suggested by the LSND and MiniBooNE experiments. For 
large Am^ effi values (>0.1 eV 2 ), the OPERA 90% C.L. upper limit on sm 2 (20 new ) based 
on a Bayesian statistical method reaches the value 7.2 x 10 -3 . 
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1 Introduction 

The OPERA experiment [1] is designed to perform an appearance search for the —> 
u T oscillations [2] in the CNGS beam [3] produced at CERN and directed towards the 
OPERA detector at the Gran Sasso Underground Laboratory (LNGS), 730 km away. A 
charged-current (CC) v T interaction in the lead-emulsion target can be identified by detect- 
ing the decay of the short-lived r lepton through particle tracking in the high-resolution 
nuclear emulsions. The observation of two v T candidate events has recently been reported 
[4-6]. The tracking capabilities of emulsions also allow to identify electrons produced in 
CC interactions of v e and therefore to search for v e appearance from v„ — > v e oscillations. 
Given the long baseline of the experiment, and the high energy of the Un, beam ((E) = 17 
GeV), OPERA has a good sensitivity for Am 2 > 0.01 eV 2 , i.e. for the LSND [7] - Mini- 
BooNE [8] allowed region (see section 4.3). Recently ICARUS, which shares the beamline 
with OPERA, severely limited this region [9]. Here we present a further constraint on these 
non-standard oscillations from the analysis of the data collected in 2008 and 2009. 

2 Detector, beam and data taking 

In OPERA, neutrinos interact in a large mass target made of lead plates interspaced with 
nuclear emulsion films acting as high accuracy tracking devices [10, 11]. This kind of target 
is historically called Emulsion Cloud Chamber (ECC). The full OPERA detector is made 
of a veto plane followed by two identical Super Modules (SM), each consisting of a target 
section and a magnetic muon spectrometer. The target sections are made of, in total, 
150,000 emulsion/lead ECC modules (or "bricks") arranged in planes, with a weight of 
about 1,250 tons, interleaved by the scintillator "Target Tracker" (TT) planes. A target 
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brick consists of 56 1-mm thick lead plates interleaved with 57 emulsion films for a total 
weight of 8.3 kg. Its thickness along the beam direction corresponds to about 10 Xo. 
Tightly packed removable doublets of emulsion films, called Changeable Sheets (CS) [12], 
are placed on the downstream face of each brick. They serve as interfaces between the TT 
planes and the bricks to facilitate the location of the neutrino interactions. 

Charged particles from a neutrino interaction in a brick cross the CS and produce 
signals in the scintillator strips of the TT. These signals are used to trigger the read-out 
and identify the brick where the interaction occurred. The brick is then extracted by 
an automated system. After development, the emulsion films are sent to the scanning 
laboratories. 

The CNGS za, beam, to which the OPERA detector is exposed, contains a small 
contamination of F„, u e , and V e . The energy spectra at the detector, as obtained from a 
Monte Carlo simulation [13], are shown in figure 1. The integrated contamination of v e and 
V e CC interactions at Gran Sasso, relative to the integrated number of CC interactions, 
is 0.88% and 0.05%, respectively. 

OPERA collected data corresponding to 17.97 x 10 19 protons on target (pot) by De- 
cember 2012 with 18941 events recorded. The analysis reported in this paper uses the data 
collected in 2008 and 2009, corresponding to 5.25 x 10 19 pot (1.73 x 10 19 and 3.52 x 10 19 pot, 
respectively) and to 5255 events recorded. The details of data taking and a comparison 
with Monte Carlo (MC) simulations for the 2008 and 2009 runs are reported in [5, 14]. 
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Figure 1. Neutrino fluxes of the different components at Gran Sasso in log scale. 
3 Emulsion scanning and search for v e interactions 

Bricks that are candidates for containing neutrino interactions are analysed following a 
complex procedure described in detail in [4, 5]. Here we just recall the main steps of the 
analysis. 
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The TT predictions are used for a large area scan of the corresponding CS films. If 
candidate tracks corresponding to the TT predictions are found in the CS, the 57 films 
contained in the brick are developed, and sent to the scanning labs. The tracks found 
in the CS are then followed upstream from film to film (scan-back) to find the neutrino 
interaction vertex. Once the vertex is found, the scanning of a volume downstream of the 
vertex (1 cm 2 in area and at least 7 films or 1.2 Xo in thickness) is performed in order to 
reconstruct all the tracks connected to the vertex and to search for decay topologies. 

The main goal is the observation of the decay of a r lepton, as a signature for a v T 
CC interaction. Moreover, the scanning also allows identifying electrons, hence v e CC 
interactions. The identification of an electron is essentially based on the detection of the 
associated electromagnetic shower. Since the size of the standard scanned volume is too 
short in the beam direction to contain the electromagnetic shower, the search for electrons 
is performed using an extended scanning volume defined by a dedicated procedure sketched 
in figure 2. All primary tracks emerging from the interaction vertex are extrapolated to 
the CS. The tracks with angles similar (A9 < 150 mrad) to that of the corresponding 
primary track (figure 2b) are searched in the CS region within 2 mm around the projected 
point. If 3 or more tracks are found in the CS, corresponding to a given primary track, an 
additional volume along the candidate track is scanned, aiming at the reconstruction of an 
electromagnetic shower (figure 2c). 



ECC 



CS 






Figure 2. Sketch for the procedure of a systematic search for v e candidates. After the reconstruc- 
tion of tracks in the standard volume (a) , all tracks emerging from the interaction vertex (the pink 
film) are extrapolated to the CS (b). If 3 or more tracks are found in the CS, corresponding to a 
given track, an additional volume along the full track length is scanned, leading to the detection of 
the electromagnetic shower (c). 



If a shower is found, the corresponding primary track becomes an electron candidate. 
The candidate track is then carefully inspected in the first two emulsion films following the 
interaction vertex. The aim is to check whether the track is due to a single particle (an 
electron) or to an e + e~ pair and so to reject electromagnetic showers initiated by the early 
conversion of a 7 from a it decay. Figure 3 shows, as an example, the reconstruction of an 
e + e~ pair from a 7 conversion in the two layers of one emulsion film. The figure illustrates 
the capability to measure tracks with micrometric resolution. To remove 7 background, 
those configurations are excluded where a track can be separated into two almost parallel 
segments more than 1 micron apart in the first or second film. 
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A significant impact parameter of the electron track with respect to the primary vertex 
would allow to identify the event as a v T CC interaction with a r — > e decay. For the present 
analysis, an upper limit of 10 /im is set on this impact parameter to select tracks originating 
from the vertex. 
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Figure 3. Side view of an e + e pair detected in an emulsion film. Note, that in a film the two 
emulsion layers are separated by a 205 microns thick plastic base. 

In addition, a scan-back procedure along the electromagnetic shower, described in [1], 
is applied. Implemented to increase the detection efficiency for a r — > e decay in the v T 
quasi-elastic interactions, it is also beneficial for the detection of v e CC interactions. 

Once the presence of an electron track is confirmed at the neutrino interaction vertex, 
the event is classified as a v e interaction. The energy of the v e candidates is estimated from 
the calorimetric measurement in the TT [14]. Moreover, a correction to the estimated en- 
ergy using a MC simulation of v e CC interactions is applied leading to the energy resolution 
for v e interactions given by the following expression: 

AE/E = 0.43 - 0.0009£ + 0.59A/VE (E in GeV). 

Among the 5255 candidate neutrino interactions collected during the 2008 and 2009 
runs, 2853 vertices were localized in the bricks, out of which 505 did not have a muon 
identified by the electronic detectors, i.e. were not classified as CC interactions. Out 
of those 505 events 19 v e candidate events were found ; 17 events were found with the 
procedure illustrated in figure 2, and the 2 remaining events were found with the scan-back 
procedure mentioned before. To illustrate the typical pattern of v e candidates, figure 4 
shows the reconstructed image of a v e candidate event, with the track segments observed 
along the showering electron track. 

The v e detection efficiency as a function of the neutrino energy is computed with a 
GEANT3 based MC simulation. The simulated events are reconstructed with the same 
algorithms as used for the data. Slight differences in the scanning strategy used along the 
years are taken into account and enter in the evaluation of the systematic uncertainty. The 
results of the simulation are shown in figure 5. The (relative) systematic uncertainty of the 
efficiency is calculated to be 10% for energies above 10 GeV and 20% for lower energies. 
Averaged over its energy spectrum, the v e beam contamination detection efficiency is £d e t = 
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Figure 4. Display of the reconstructed emulsion tracks of one of the v e candidate events. The 
reconstructed neutrino energy is 32.5 GeV. Two tracks are observed at the neutrino interaction 
vertex. One of the two generates an electromagnetic shower and is identified as an electron. In 
addition, two showers from 7 conversions are observed (overlapping in this projection), starting 
from 2 and 3 films downstream of the vertex. 



(53±5)%. For neutrino energies smaller than 30 GeV (20 GeV) it is e det = (43±5)% ((35± 
4)%). 
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Figure 5. Detection efficiency of v e events as a function of the neutrino energy, obtained from MC 
simulations. The error bars show the estimated systematic uncertainties. 



Apart from the v e background associated to the beam contamination (see section 4), 
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two main sources of background are considered for the v e search: (a) 7r° misidentified as 
electron in neutrino interactions without a reconstructed muon; (b) v T CC interactions 
with the decay of the r into an electron. 

Background (a) occurs if an e + e~ pair appears to be connected to the interaction 
vertex and cannot be distinguished from a single particle in the first two emulsion films 
after the vertex or if one branch of the pair has a very low energy and remains undetected. 
This background was evaluated directly from the data. We have analysed the 7 conversion 
patterns observed in a sample of 1106 neutrino events and studied the probability that 
they could be selected as electron candidates. The effect of ir° decays to a Dalitz pair is 
estimated to be negligible. By normalizing to our sample, the estimated background (a) is 
0.2 ± 0.2 events. 

Background (b) was computed by MC simulation assuming the three-flavour — > v T 
oscillation at maximal mixing and Am 2 = 2.32 x 10~ 3 eV 2 [15]. Requesting that the r decay 
topology was not detectable (no kink, no significant impact parameter of the electron track 
with the interaction vertex), we calculated this background in our sample to be 0.3 ±0.1 
events. 

The total amount of the considered background for the v e CC interaction search is 
0.4 ± 0.2 events. 

4 Oscillation analysis 

4.1 Background to — > v e appearance 

A signal for Un — > v e oscillations should appear as a significant excess of electron events 
with respect to the expected background, mainly due to v e and V e CC interactions from 
the beam contamination. A detailed evaluation of this number was performed, starting 
from the fluxes of the different beam components presented in figure 1. The systematic 
uncertainty on the v e flux was estimated to be 4% on its normalization and 3-4% on its 
energy dependence [16]. We conservatively applied a 10% systematic uncertainty on the 
overall v e beam contamination. The fluxes were weighted with the CC cross sections and 
the energy dependent detection efficiency. An additional inefficiency of 6% (for 2008) and 
3% (for 2009) was introduced to reflect the lower film quality of a small fraction of the 
bricks. Taking into account the target mass and the pot corresponding to our data, we 
expect to observe 19.4 ± 2.8 (syst) u e events from the beam contamination in the full 
energy range. Together with the backgrounds (a) and (b) discussed above, we expect 
19.8 ± 2.8 (syst) background v e events. This number is in agreement with the 19 observed 
candidate v e events and therefore the room for oscillations is reduced. In the following we 
analyse two scenarios. 

4.2 Three- flavour mixing scenario 

A non-zero #13 has recently been reported by several experiments [17-20]. Using the 
following oscillation parameters [15] : sin 2 (2#i3) = 0.098, sin 2 (26>23) = 1, Ara| = Arn^j = 
2.32 x 10 -3 eV 2 , also assuming 5qp — and neglecting matter effects, 1.4 oscillated v e CC 
events are expected to be detected in the whole energy range. 
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Energy cut 


20 GcV 


30 GeV 


No cut 


BG common to 
both analyses 


BG (a) from vr° 


0.2 


0.2 


0.2 


BG (b) from r —> e 


0.2 


0.3 


0.3 


v e beam contamination 


4.2 


7.7 


19.4 


Total expected BG in 3-flavour oscillation analysis 


4.6 


8.2 


19.8 


BG to non-standard 
oscillation analysis only 


v P via 3-flavour oscillation 


1.0 


1.3 


1.4 


Total expected BG in non-standard oscillation analysis 


5.6 


9.4 


21.3 


Data 


4 


6 


19 



Table 1. Expected and observed number of events for the different energy cuts. 



Figure 6 shows the reconstructed energy distribution of the 19 v e candidates, compared 
with the expected reconstructed energy spectra from the v e beam contamination, the os- 
cillated v e from the three- flavour oscillation and the background (a) and (b), normalized 
to the pot analysed for this paper. To increase the signal to background ratio a cut E < 20 
GeV is applied on the reconstructed energy of the event, which provides the best figure of 
merit on the sensitivity to #13. Within this cut, 4.2 events from v e beam contamination 
and 0.4 events from the backgrounds (a) and (b) are expected, while 4 events are observed. 
The numbers are summarized in table 1. The number of observed events is compatible 
with the non-oscillation hypothesis and an upper limit sin 2 (2#i3)< 0.44 is derived at the 
90% Confidence Level (C.L.). 
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Figure 6. Distribution of the reconstructed energy of the v e events, and the expected spectrum 
from the different sources in a stack histogram, normalized to the number of pot analysed for this 
paper. 
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4.3 Non-standard oscillations 

Beyond the three- neutrino paradigm, some possible hints for non-standard effects have 
been reported, in particular by the LSND [7] and MiniBooNE [8] experiments. We have 
used OPERA data to set an upper limit on non-standard Un — > v e oscillations. 

We used the conventional approach of expressing the za, — > v e oscillation probability 
in the one mass scale dominance approximation, given by the following formula with new 
oscillation parameters 9 new and Am^ ew : 

P u ^ Ue = sm 2 (29 new ) ■ sin 2 (1.27 Am 2 new L(km)/E(GeV)) 

Note however that this approach does not allow a direct comparison between experiments 
working in different L/E regimes [21]. 

The flux at the detector, normalized to the integrated statistics used in our anal- 
ysis, is weighted by the oscillation probability, by the CC cross-section and by the energy 
dependent detection efficiency, to obtain the number of v e CC events expected from this 
oscillation. 

As the energy spectrum of the oscillated v e with large Am^ etll (>0.1 eV 2 ) follows the 
spectrum of z^, which is basically vanishing above 40 GeV (see figure 1), a cut on the 
reconstructed energy is introduced. The optimal cut on the reconstructed energy in terms 
of sensitivity is found to be 30 GeV. We observe 6 events below 30 GeV (69% of the 
oscillation signal at large Am 2 eiu is estimated to remain in this region), while the expected 
number of events from background is estimated to be 9.4 ± 1.3 (syst) (see table 1). Note 
that we choose to include the three-flavour oscillation induced events into the background. 
In this case, the oscillation probability does not contain the #13 driven term. 

The 90% C.L. upper limit on s'm 2 (26 new ) is then computed by comparing the expec- 
tation from oscillation plus backgrounds, with the observed number of events. Since we 
observed a smaller number of events than the expected background, we provide both, the 
Feldman and Cousins (F&C) confidence intervals [22] and the Bayesian bounds, setting a 
prior to zero in the unphysical region and to a constant in the physical region [23]. Un- 
certainties of the background were incorporated using prescriptions provided in [15]. The 
results obtained from the two methods for the different C.L. are reported in table 2. We 
also quote our sensitivity calculated assuming 9 observed events (integer number closest to 
the expected background). 

Given the underfluctuation of the data, the curve with the Bayesian upper limit was 
chosen for the exclusion plot shown in figure 7. For convenience, results from the other 
experiments, working at different L/E regimes, are also reported in this figure. For large 
Am 2 ew , values the OPERA 90% upper limit on sin 2 (2# ne „,) reaches the value 7.2 x 10~ 3 , 
while the sensitivity corresponding to the pot used for this analysis is 10.4 x 10 -3 . 

5 Conclusions and perspectives 

First results of a search for — > v e oscillations with the OPERA experiment at the Gran 
Sasso Underground Laboratory have been presented. The experiment searched for the 
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Upper limit 


Sensitivity 




C.L. 


r &C 


Bayes 


r &C 


Bayes 


Number of oscillated 


nrW 

yU/o 


3.1 


4.5 


6.1 


6.5 


v e events 


95% 


4.3 


5.7 


7.8 


7.9 




99% 


6.7 


8.2 


10.7 


10.9 


sm 2 (29 new ) at 


90% 


5.0xl(T 3 


7.2xl(T 3 


9.7xl0~ 3 


10.4xl0~ 3 


large Am 2 


95% 


6.9xl0~ 3 


9.1xl0~ 3 


12.4xl0~ 3 


12.7xl0~ 3 




99% 


10.6xl(T 3 


13.1xl0~ 3 


17.1xl0~ 3 


17.4xl0~ 3 



Table 2. Upper limits on the number of oscillated v e CC events and sin (20 netu ), obtained by the 
F&C and Bayesian methods, for C.L. 90%, 95%, 99%. The sensitivity is computed assuming that 
the number of observed events is 9, which is the closest integer to the 9.4 expected background 
events. 




Figure 7. The exclusion plot for the parameters of the non-standard — > v e oscillation, ob- 
tained from this analysis using the Bayesian method, is shown. The other limits shown, mostly 
using frequentist methods, are from KARMEN (p^ — > v e [24]), BUGEY (V e disappearance [25]), 
CHOOZ (V e disappearance [26]), NOMAD -> v e [27]) and ICARUS -> v e [9]). The regions 
corresponding to the positive indications reported by LSND (77 M — > v e [7]) and MiniBooNE (v^ — > 
v e and v ^ — > v e [8]) are also shown. 

appearance of v e in the CNGS neutrino beam using the data collected in 2008 and 2009, 
corresponding to an integrated intensity of 5.25 x 10 19 pot. The observation of 19 v e 
candidate events is compatible with the non-oscillation expectation of 19.8±2.8 events. 

The current result on the search for the three-flavour neutrino oscillation yields an 
upper limit sin 2 (26>i 3 ) < 0.44 (90% C.L.). 

OPERA limits the parameter space available for a non-standard v e appearance sug- 
gested by the results of the LSND and MiniBooNE experiments. It further constrains the 
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still allowed region around Am^ M = 5 x 10 -2 eV 2 . For large Am 2 eu) values, the 90% C.L. 
upper limit on sin 2 (20 ne iu) reaches 7.2 x 10 -3 . This result is still affected by the statistical 
underfluctuation, the sensitivity corresponding to the analysed statistics being 10.4 x 10 -3 . 
A Bayesian statistical treatment has therefore been adopted for determining the upper 
limit. 

Various improvements are expected for the future. The statistics will be increased 
by a factor of 3.4 by completing the analysis of the collected data. The reconstructed 
energy resolution will be improved when the calorimetric measurement in the TT will be 
complemented by following the hadron tracks and the electron showers in the downstream 
bricks. 

With the increase in sample size and the improvements in the analysis, the effect of a 
possible statistical underfluctuation of the background will be reduced and OPERA should 
then be able to access the parameter region comparable to its sensitivity below sin^(2^ netL) ) 
= 5.0xl0~ 3 . 
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